www.rspsciencehub.com 





BHA WS’ VO, 


INTERNATIONAL RESEARCH JOURNAL ON 
ADVANCED SCIENCE HUB 


RSP SCIENCE HUB 


(The Hub of Research Ideas) 
Available online at www.rspsciencehub.com 


Volume 03 Issue 04 April 2021 


e-ISSN : 2582 - 4376 
Open Access 


Stress Intensity factors of a semi-elliptical crack In a Compressor blade of 


dovetail joint through numerical approach 


Dr.Chandramouli Harige Ramaswamy’, Dr. Chakradhar Goud * 
‘Department of Mechanical Engineering, Shri JJT University, Jhunjhunu, Rajasthan, India. 
’Professor & Principal, Springfields Engineering College, Hyderabad, Telangana. India. 


harigechandramouli@ gmail. com’, cgsakki@ gmail. com 


Abstract 


Induced fatigue has led to a large number of deviations of service failures in a wide range of applications, 
where the deterioration of properties occurs due to the slipping between surfaces at a small displacement. 
In this article, a numerical simulation of a semi-elliptical crack in the dovetail joint of the compressor was 
performed to understand the destruction of the blade due to the stress contour in the contact area. A 
three-dimensional crack is modelled using the ANSYS Finite Element Analysis program, and the blade 
profile is modelled by CATIA V5R20 software. The study considers a model of an open semi-elliptical 
crack, and the parameters of the movement of the crack tip are estimated using three-dimensional 
singular elements of the crack tip in ANSYS. Initially, a three-dimensional finite element model of a 
compressor blade with a dovetail contact area was simulated in CATIA V5 and analysed in the ANSYS 
working environment to determine the maximum principal stresses in the crack-free region. After the 
crack front has been precisely identified and confirmed, fracture modelling is used to determine the 
impact of crack size with one side, as well as rotation speed, on fracture parameters (stress intensity 


factors). 
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1. Introduction: 

The root of the dovetail of the compressor blades 
is the case when the induced fatigue failure is a 
fairly common occurrence. Under the influence of 
centrifugal loading and vibration, fatigue occurs 
between the root of the blade, the disk and the 
stresses in the vicinity of the contact, which must 
be accurately estimated to predict the fatigue life. 
Stress assessment is a big problem for designers 
because of the complex geometry, load, and non- 
linear _behaviour.[1-4].Although — considerable 
attention has been given in recent years to the 
analysis of failures of turbine blades with cracks, 
and also due to advances in 3D fracture analysis, 
efforts have been made to characterize the 
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propagation of cracks in crack blades. In this 
study, the stress intensity factor is analysed based 
on the location where the theory of maximum 
principal stress 1s considered, and semi-elliptical 
cracks are modelled using the finite element 
method (FEM).[5-8]. P. Papanikos et al. conducted 
a three-dimensional nonlinear analysis by the finite 
element method of this area of the dovetail in the 
disk assemblies of the aircraft engine compressor 
using contact elements. Lucian witek (conducted a 
failure analysis of an aircraft engine turbine disk. 
Nonlinear finite element analysis was used to 
determine the stress state of the disk / blade 
assembly under various rotation conditions.[8-10]. 
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Crack initiation sites 
(b) 


Fig.1. Crack initiation sites on the compressor blade 


2. Description of the problem 

In this study, a semi-elliptical crack is modelled 
based on maximum principal stress and Von mises 
stresses region on the compressor blade profile 
without crack. Then maximum Stress intensity 
factors for semi-elliptical crack at various contours 
are calculated for different crack lengths and 
rotation of the disc. Titanium alloy blade material 
is considered for analysis which is a good 
corrosion resistant, high chemical stability and 
excellent mechanical properties. 

3. Numerical simulation of compressor blade in 


13.11 


Blade Disc A 


dovetail joint 

3.1. Geometric modelling 

The geometric modelling was carried out using 
CATIA-V5. In view of disc symmetry only one 
sector of the blade was modelled for analysis as 
shown in the Fig. 2 and Fig.3. The element is 
defined by 10 nodes having three degrees of 
freedom at each node: translations in the nodal x, 
y, and z directions. For modelling semi elliptical 
crack, major radius of 0.05 mm and minor radius 
of 0.04 mm are considered. 


13.05 


Blade B 


All dimensions in mm 





Fig.2.Geometric details of a compressor Blade dovetail joint 


International Research Journal on Advanced Science Hub (IRJASH) 


62 


www.rspsciencehub.com 


3.2. Finite Element modelling 

ANSYS software was used for the Discritization of 
the 3D model. Using SOLID187 components, the 
blade and disc were discretized. SOLID187 is a 
10-node tetrahedral 3D structural element used to 
3.3. Contact properties 

To perform non-linear analysis, a symmetric 
contact was established from a solid to a solid 
contact of the blade and disk. CONTAI74 and 
TARGET170 were used to represent the contact 
and target surfaces, respectively. The surface of 
the blade was modelled as a contact, and the 
surface of the disk was modelled as a target 
surface, as shown in Fig. 4. CONTAI174 is used 
for three-dimensional structural and pair contact 
analysis. This element can be located on the 
surfaces of three-dimensional solid or shell 
elements that do not have any middle lateral nodes. 





Fig. 3. Geometric details of a compressor blade 
dovetail joint 3D modeling 


3.4. Loads and Boundary Conditions 

To run the analysis only centrifugal load was 
considered. The boundary constraints are one of 
the most important factors in Finite element 
analysis. The boundary condition for compressor 
blade with dovetail joint is as shown in the Fig.6. 
The rotor and blade have cyclic symmetry about 
the axis hence only one sector of the model is 
considered for analysis. The bottom portion of the 
disc was constrained radially and the side faces of 
the disc were constrained using symmetry 
boundary conditions. 
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model solid structures in three dimensions. Figure 
4 depicts the finite element modelling of a 
compressor blade dovetail joint assembly and a 
semi-elliptical crack. 


When the surface of an element penetrates one of 
the target elements of a segment on a given target 
surface, contact occurs. TARGET170 is a three- 
dimensional target surface for linked contact 
elements that is used to describe various three- 
dimensional target surfaces. The target surface is 
discretised by a set of elements of the target 
segment and connected to the corresponding 
contact surface through a common real constant 
set. In Fig.5, it’s shown that there are total 113377 
nodes, number of contact elements are 12113, 
number of solid elements are 66942 and number of 
total elements are 79056. 





(b) 


Fig.4. Finite element model (a) compressor 
dovetail joint (b) Semi-elliptical crack 


4. Results and discussions 
In this section, the results obtained from carrying 
out non-linear static analysis on 3-Dimensional 
dovetail joint of compressor are presented. 
4.1Analysis of Compressor Blade dovetail joint 
without Crack 

In Fig.7. and Fig.8.show both the von-mises 
Stress and the maximum principal stress at the 
contact region, which gives the identification of 
crack location at the contact region of the dovetail 
joint. 
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A-Coupling constraint 
B-Centrifugal force 
C-Fixed constraint 
D-Radial Constraint 


Fig.5. Modelling ofDovetail jointContact 
assembly 


924, 76 Max 
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116.62 

58,315 
0.009248 1 Min 


Fig.7. Maximum Von-mises stressat the 
contact region without crack 


4.2 Analysis of Compressor Blade dovetail joint 
with Crack 

Based on the results of the 3D analysis, high 
stressed location near the blade and disc contact 
region which is the probable location of initiation 
of crack was identified. In order to predict the 
stress Intensity factor, a semi-elliptical crack was 
introduced in the region where maximum von- 
mises stress and maximum principal stress occurs, 
which is shown in the fig.9 and Fig.10. The 
KCALC command used in plain strain option to 
calculate the maximum stress Intensity factor (K1) 
at the crack region, which is predominant, 1s 
shown in the Fig.11. From Fig.9, it’s observed that 
von-mises stress is very high at the crack tip (823 
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Fig. 6. Load and boundary Conditions of dovetail 


assembly 


543.81 Max 
Af5.O7 
406,14 


340.3 

ef2.47 
24.03 

136.8 

68,964 
1.1287 

-66. 706 Min 





Fig.8. Maximum Principal Stress at the contact 
region without crack 


MPa) because of singular elements and 
concentrated plastic zone at the crack tip. 
Maximum principal stress of 1339 MPa is 
observed from the fig.10, which is higher than the 
ultimate stress value due to the crack growth rate. 
Fig. Illrepresents stress Intensity factors K1 for 
mode I. 

4.3 Evaluation of stress Intensity factor (SIF) 
for different modes in compressor disc. 

The results obtained from carrying out non-linear 
static analysis on 3D dovetail joint of compressor 
are presented here. From Fig.12, the results 
obtained are presented in terms of modes of I, II 
and III stress intensity factors (K1,K2 and K3) 
variations over various crack depth for semi- 
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elliptical for a compressor speed of 8000 r.p.m (@ 
=837.75 rad/s) for six contours. In Fig. 12. the 
values of stress Intensity factors for K1, K2 and 


623.05 Max 
#31,6 

640,15 

345, / 

457.25 


363.0 

274,36 

1382.91 

91.456 
0.0090589 Min 


Fig.9 Maximum Von-mises stress at the crack 
region 
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30.302 
49.348 
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Fig.11 Maximum stress Intensity factor(K1) at 
the crack tip region 

Conclusions 

An effort has been made in the present study to 
understand the analysis of semi-elliptical crack in a 
compressor blade dovetail joint assembly through 
numerical simulation. Three-dimensional non- 
linear finite element of a compressor blade 
dovetail joint assembly is carried out to evaluate 
the maximum stress in the contact region. The 
semi-elliptical crack analysis was carried out to 
find out the maximum stress Intensity factors for 
three fracture modes of different crack lengths. 
The SIF increases as the length of the crack depth 
increases. The Stress Intensity factor for a given 


Stress Intensity factor(MPa-mm”"0.5) 
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K3 are plotted for fracture mode 1, mode 2 and 
mode 3 respectively. 


1339.3 Max 
1183.1 
1026.9 
8/0.67 
714,48 
558.28 
402.08 
245.89 

89.69 
-66.50 7 Min 





Fig.10 Maximum Principal Stress at the crack 
region 
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Fig.12 Values of K1, K2 and K3 for six contours 


compressor speed is evaluated through the induced 

semi-elliptical crack 1s to predict the crack growth 
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